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The Effect of Primordial Non-Gaussianities on the Seeds of Super-Massive Black Holes 
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The origin of the seeds which develop into the observed super-massive black holes at high redshifts 
may be hard to interpret in the context of the standard A CDM of early universe cosmology based 
on Gaussian primordial perturbations. Here we consider the modification of the halo mass function 
obtained by introducing skewness and kurtosis of the primordial fluctuations. We show that such 
primordial non-Gaussianities constrained by the current observational bounds on the nonlinearity 
parameters of f nl and grjvL are not effective at greatly increasing the number density of seeds which 
could develop into super-massive black holes at high redshifts. This is to be contrasted with the role 
which cosmic string loops could play in seeding super-massive black holes. 

PACS numbers: 98.80.Cq 


I. INTRODUCTION 

Super-massive black holes (SMBHs) with masses ex¬ 
ceeding 10 6 M 0 (M q denotes the solar mass) have recently 
attracted a lot of attention [IH331- It is now believed that 
each galaxy contains at least one super-massive black hole 
which forms from the accretion of gas about massive seed 
objects. The origin of the seeds which cause the forma¬ 
tion of SMBHs is still somewhat of a mystery [HI IT5] . 
According to the standard paradigm of early universe in 
which the primordial cosmological fluctuations are ap¬ 
proximately Gaussian and have an almost scale invariant 
spectrum, nonlinearities form only at fairly late times and 
there may not be enough time to produce the nonlinear 
massive seeds which seed SMBHs of mass greater or equal 
to 10 9 M 0 at redshifts of 6 or higher (of which more than 
40 candidates are now known 0E|). There are three 
types of candidate seeds for SMBHs, namely Population 
III stars with masses in the range between 10 2 Mq and 
10 3 Mq, dense matter clouds with masses between 10 3 A4q 
and 10 6 Mq, and compact objects of mass between 10 2 A4q 
and 10 4 Mq formed by the collision of old stellar clusters. 

Recently m it has been shown that cosmic string loops 
which result from a scaling solution of strings formed dur¬ 
ing a phase transition in the very early universe lead 
to an additional source of compact seeds. The number 
density of string-induced seeds dominates at high red¬ 
shifts and can help trigger the formation of the observed 
super-massive black holes. Cosmic string loops form a 
special type of non-Gaussian density fluctuation field. 
Non-Gaussianities in the primordial density perturbation 
field are an inevitable consequence of the nonlinearities of 
the Einstein field equations. However, in minimal models 
of matter and in the context of an inflationary origin of 
the density fluctuations such non-Gaussianities have a 
very small amplitude. Larger amplitudes can be obtained 
in non-minimal models of inflation nz] and in some al- 
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ternatives to inflation such as the “Ekpyrotic” scenario 
[T81J and in the “matter bounce” [T5]. In such models, 
the non-Gaussianities are usually parametrized in terms 
of the skewness and kurtosis which characterize the de¬ 
viations from Gaussianity in the three and four point 
functions. Skewness and kurtosis are a good character¬ 
istic of non-Gaussianities emerging from a distortion of 
originally Gaussian fluctuations. In this paper we ask 
whether such non-Gaussianities can play a similar role as 
cosmic string loops in the triggering of the formation of 
nonlinear seeds for SMBHs at high redshifts. 

Primordial perturbations which originate as quantum 
fluctuations of a scalar field and which seed cosmic mi¬ 
crowave background (CMB) fluctuations and structure 
formation of the Universe lead to non-Gaussianities which 
are well characterized by skewness and kurtosis (in con¬ 
trast to string loop-induced fluctuations which are intrin¬ 
sically non-Gaussian and hence not well characterized by 
the three and four point functions only. Nevertheless, the 
conventional non-Gaussianities resulting from a deforma¬ 
tion of an initially Gaussian process can yield information 
about the early stages of the evolution of our Universe. In 
particular, they can allow discrimination between various 
models of inflation and their alternatives. 

In general, the three and four point functions are char¬ 
acterized by an amplitude and a shape. Here, we will 
focus on local type non-Gaussian perturbations which can 
be written in the following form [201 GQ] 

C = Cg + g/jVL(Cc - (Cg)) + 5 9nl(g■ (!) 

where £ is the primordial curvature fluctuation variable, 
£<3 is the Gaussian part and /nl and gwL are nonlinearity 
parameters. These nonlinearity parameters /a tl and gNL 
lead to non zero values of the skewness and kurtosis 
of the primordial non-Gaussianities in the Probability 
Density Function (PDF), respectively IHHH] , Current 
observations constrain the magnitude of the parameters 
f nl and g NL [26] to be = -16 ± 79 and g'^f 1 = 
(—9± 7.7) x 10 4 (both at 68% confidence level). Here, the 
superscripts stand for the shape, “equilateral” in the first 
case and “local” in the second. 
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In this paper we take into acount the modification of the 
halo mass function in the presence of skewness and kur- 
tosis of local type (following the method of j22Jj27|) and 
determine the effects of these primordial non-Gaussianities 
on the mass of the dark matter halos at different redshifts 
(up to ^ = 20) compared with the situation in a purely 
Gaussian model. We demonstrate that the constraints 
which come from the current observational data on the 
magnitude of the nonlinearity parameters are already 
strong enough to prevent the non-Gaussianities from elim¬ 
inating the difficulties that the standard A CDM of early 
universe cosmology may face in terms of explaining the 
origin of the observed super-massive black holes at high 
redshifts. 

The paper is organized as follows. First we start with a 
short review of the halo mass function in the A CDM 
model. In Section III we study the modifications of 
the halo mass function resulting from introducing non¬ 
vanishing skewness and kurtosis of the primordial density 
perturbations. We evaluate the result for skewness and 
kurtosis of equilateral shape. We show that this type 
of primordial non-Gaussianities obeying the current ob¬ 
servational bounds on the nonlinearity parameters f^L 
and gpjL is not effective at greatly enhancing the number 
density of nonlinear seeds for SMBHs at high redshifts. 
We conclude with a discussion section. 


II. HALO MASS FUNCTION IN THE A CDM 
MODEL 

The linear relative matter density fluctuation can be 
written in terms of the primordial curvature fluctuations 
Q(k) on uniform energy density hypersurfaces by [22] [25] 

p ;.2 

S(k,z) = —^T{k)D{z)ak) , (2) 

where k denotes comoving wavenumber, z is the cosmo¬ 
logical redshift, fl m0 is the present density parameter, H 0 
is the Hubble constant, D(z ) is a linear growth function 
and T(k) is a transfer function. Hence, the linear matter 
power spectrum Ps(k) is given by the curvature power 
spectrum Pq via 

(S(k 1 ,z) 6 (k 2 ,z)) = (2ir 3 )P s (ki,z)5 3 (k 1 +k 2 ) 

Ps(ki,z) = A(ku z ) 2 Pf(fci) , (3) 

where 

9k 2 

A(k,z) = ——iT{k)D{z). (4) 

01 LuiO-LIq 

The smothed density fluctuation on a given length, R, 
is defined by 


/ r1 3 k 

—W R (k) 6 (k,z), (5) 


with W R (k ) being the Fourier transform of a window 
function 


W R (k) 


3 (sin kR — kR cos kR) 
(kR ) 3 


( 6 ) 


The variance in mass on a momentum scale k with a 
top-hat filter with a radius that encloses the mass is (for 
an infinite total spatial volume) given by 


a R = j ~j 7 W R(k) A(k,z) 2 V c (k) 

= I fw 2 R (k) V s (k) , ( 7 ) 


where the dimensionless power spectra V are related to 
the dimensional ones P via 


■pm = ( 8 ) 

The halo mass function based on the Press-Schechter 
theory [29] is given by [30] 


= 


-dM^ d 


= -dM 


ar — / dvF(v) (9) 
M dM J Sc/aR 

2 p r vl dlog a R 


and 


dvF( v) 




( 10 ) 


where p is the mean mass density of the universe (back¬ 
ground density) and v c = The number 5 C is the 
threshold for collapse, and we will use the value S c = 1.86 
which corresponds to neglecting the effects of dark energy. 

According to linear perturbation theory, the density 
contrast a R grows linearly in the scale factor before the 
contribution of dark energy to the equation of state of 
matter becomes important. The above formulas then 
show that for a Gaussian spectrum of fluctuations, the 
number density of nonlinear seeds falls off exponentially at 
high redshifts, thus making it difficult to account for the 
origin of the nonlinear seeds which are needed to explain 
the origin of SMBH. We will now investigate whether 
the addition of skewness and kurtosis can improve the 
prospects for high redshift SMBH formation. 


III. MODIFICATION OF THE HALO MASS 
FUNCTION BY PRIMORDIAL 
NON-GAUSSIANITIES 


In this section we focus on the local type non- 
Gaussianities and study the modifications of the halo 
mass function when introducing primordial skewness and 
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kurtosis (see also J3T] [32] for other studies of the effects of 
non-Gaussianities on structure formation). We take the 
primordial curvature fluctuations to be given by 0 in 
terms of a Gaussian distribution. Based on this relation 
the two, three and four point functions become 

(C(fci)C(fc 2 )) = (27r) 3 P c (fc 1 )<5 (3) (fc 1 + fc 2 ) (11) 
(C(fci)C(fc 2 )C(fc 3 )) = (2tt ) 3 ®/ JVi (P < (fc 1 )P c (fc 2 ) 

+2 perms.) S^ 3 \k i + k 2 + fc 3 ) 
_ _ _ _ 54 

<C(* 1 )C(te)C(*3)C(*4)) = (2ir) 3 -g N L(Pdki) 

P c (/c 2 )P c (fc 3 ) + 3 perms.) 

6 ^ (fci + fc 2 + k 3 + hi) . 

In order to consider the effects of primordial non- 
Gaussianities on the smoothed density fluctuations, we 
can define the n-th central moment of the Probability Den¬ 
sity Function (PDF) P(<5 r)cMr in the following standard 
way 

/ OO 

6 n R F(S R )d 6 R . (12) 

-OO 

In addition, we can define the reduced p-th cumulant S p 
in this form 


where 


P } (^)r 1 ’ 


(Sr) c = 0 , {S 2 r ) c = cr 2 R , etc. 


(13) 


(14) 


with 


c 3 = - S 3 (R)a R , , c 4 = S 4 (P)4 , .... (19) 


This technique allows us to calculate the non-Gaussian 
PDF of the density field in terms of the non zero values of 
the skewness and kurtosis of the primordial fluctuations. 
Inserting (16) into (15) and making use of (19) we obtain 


F{v)dv = < ^— exp(—i' 2 /2)[l + h 3 {v) (20) 


V 2 n 
I.S 3 (R)<t r 
2 


+ ^^^) a H e (y) + U^^) 3 H 9 (v) 


6 


6 

1 S 3 (R)a Ry 

6 


6 


Sa(R)(T 2 r 1 SdR)v 2 Rs2TT / A 

+ 24 24 } H&{V) 

l 1 (S 4 (R)a R3 

+ 6 ( 24 } H ^>+ -] > 

where we are neglecting terms of higher order in Sn(R) 
(from 00), © and 05])), and where 


S 3 (R) = 


6 % [ ^WR+h) AikdVdh) 


5 a 


r J k\ 

' dk 2 


( 21 ) 


x J A(k 2 )V c (k 2 ) 


rl dp 


12 


'-1 


W R (k 12 ) A(k 12 ) (l + , 


and 


-Gaussian PDF of matter 

SdR) 54 

9nl / 

12), we can build up the 

b * [K) 25 

4 J 


-r^Wnih) AikdVdh) 

K\ 


( 22 ) 


PDF from the cumulants and the Gaussian distribution 
by using the Edgeworth expansion. This technique gives 
the non-Gaussian PDF of the density field in terms of 
derivatives of the Gaussian PDF and reduced cumulants 


x j ^ W R (k 2 ) A{k 2 )Vdk2) 
x j ^t 3 W *( k 3) Ak 3 )V c (k 3 ) 


F(v)dv = Fq(v) dv + Y C ^F < £ n \v) du 


m =3 


rl dp i 2 f l dpi 23 


(15) 


d(f> 13 W R (k\ 23 ) A{k\ 23 ) 


with 


x 1 + 


2 7-i 2 

Pdki2 3 ) P((ki 23 ) Pf(fc 123 ) 


Pdk 1 


Pdk 2 


Pdks) 


Fg(v) = 


V 2 d 


exp{—v 2 / 2) 


(16) 


4 m) M = = (-lr^MFeM , 

where Fed) is the Gaussian PDF and H m (v) are the 
Hermite polynomials 

Hi(v) = v, H 2 {v) — v 2 - 1 , H 3 (y) = d - 3is 
H 4 d) = v 4 - 6i/ 2 + 3 , . (17) 

The coefficients c m are given by 


/ OO 

H m (v)F(y)dv , 

-OO 


(18) 


and where 

f^ij — COS Oij , 

fc 12 = \Jk 2 + k 2 + 2 k\k 2 pi 2 and 
k\ 23 = k 2 + k 2 + k 2 + 2 k\k 2 pi 2 + 2 k\k 3 p\ 3 

+2k 2 k 3 ( J (1 - p\ 2 ){l - M 13 ) cos 013 + Mi 2 Mi3 


(23) 


1/2 


Thus, the modified halo mass function based on the 
Press-Schechter formula with a non-Gaussian PDF of the 
smoothed density field is given by (from Eqs. 0 and 
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( 20 )) 


m (M,^M = -dMgA| 0 ° dvF[y) 


= -dM 


2 P r C 

j« exp[ -yi 


6 c/< 7 R 

d log era 
dM 


k-space filter which is particularly useful for theoretical 
arguments and is the equivalent to the top-hat filter in 
(24) Fourier space [28] , 

= (29) 




1 „ , > (S 3 <trY 1 




+ 2^8(^c) 


m 

2 \ 3 n 


K-i-) 

+\n,M + \hm ( ,S *'('"') 


dM \ 6 / 


dM V 24 / 


, 1 ^ d f S 3 a R \ 2 1 ^ d y S^crf^ 

+ o^5(^c)wT 7 —+ o^7(^c)yT7 <Ti 


dM y 6 y 
, 1 „ / / S 3 a R \ i 

+ 6 Hs ^ c dM V 6 J 

2 , d (S,al V 

+ 6^ llK) dM 


dM V 24 y 


here v c = and d c is the threshold for collapse. 

In what fellows we consider two important forms of 
non-Gaussianity for cosmological observations, namely 
the local form and the equilateral form. The local form of 
the bispectrum requires that one of the three momentum 
modes exits the Hubble radius (e.g. in the context of an 
inflationary universe) much earlier than the other two, i.e 
fci <C &2 — k 3 . Taking this form of the non-Gaussianity 
for skewness and kurtosis we get 

S 3 (R) = f A(h)Vdki) ,(25) 

5 cr R J ki 

and 


In real space this takes the form 

. . 3 ,x._o, . Ixl Ixl Id, 

Wr(x) = — |-| 3 (sin— - — cosy) , (30) 

where the volume is taken to be V = 6ir 2 R 3 . Therefore, 
we obtain 

1 O /- Q o 

S 3 (R)a R ~ —fNLy/Vdl/R) ^ — x lO~ 3 f NL (31) 
and 


Sa{R)cj 2 r ~ —9nlP(; 0-/R) — ~^r x 10 ‘g NL 


(32) 


where we used the approximations n s = 1 and A s = 
2.4 x 10~ 9 . This analytical analysis shows that only 
primordial non Gaussianities with magnitudes /nl > 10 2 
and gNL > 10 6 yield an important contribution to the 
modified halo mass function of Eq. (24). 


Substituting Eq. (271 into the modified halo mass func¬ 


tion Eq. (24) we can compute the mass of nonlinear halos 


with a mean separation d as a function of redshift. Having 
such haloes is a necessary condition for the formation of 
super-massive black holes. Since there is evidence that ev¬ 
ery galaxy harbors a SMBH, we are interested in the value 
d = dg a i corresponding to the mean comoving separation 
of current galaxies 




(33) 


S *(R) = Aki)Vdki) 

x J 2 ) A(k 2 )V c (k 2 ) . (26) 

For the equilateral form of non-Gaussianity, i.e k 3 = 
k 2 = k 3 = ki = fc, we have 

S 3 (R) = y^r/ A(k)V c (k) (27) 

and 

S4(R) = fw R (k) A(k)Vdk)) 2 • (28) 

We can now estimate the magnitude of the skewness and 
kurtosis in the equilateral limit. Let us consider a sharp 


By using Eq. (33) and the CAMB code |33j we can 
plot the diagram of dark matter mass of halos with mean 
separation d ga i as a function of redshift for different values 
of the non-Gaussianity parameters /nl and gNL and 
assuming equilateral shape. The results are shown in 
Figure 1 (see also [31} for an analysis of the halo mass 
function for Gaussian fluctuations). 

Our results show that the effect of primordial non- 
Gaussianities given by f nl and gNL with values of these 
parameters consistent with the current observational 
bounds cannot change the number distribution of high red¬ 
shift nonlinear dark matter halos in a way which will have 
an important impact on the formation of super-massive 
black holes at high redsliifts. To have an important ef¬ 
fect at a redshift z = 20 values of the non-Gaussianity 
parameters larger than /nl > 10 2 or gNL > 10 6 would 
be required. 
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FIG. 1: Nonlinear dark matter halo mass (vertical axis, in solar mass units) as a function of redshift (horizontal axis, up to 
z = 20) obtained by using the CAMB {33] code for different values of the non-Gaussianity parameters Jnl and gNL and assuming 
equilateral shape. The mass shown is the mass which is turning nonlinear on a co-moving length scales which corresponds to the 
separation d ga i = IMpc of galaxies. We have made use of the value = 10 23 Mq. We have shown curves corresponding to the 
lowest values of /j ml and gNL for which the non-Gaussianities have an important effect. 


IV. CONCLUSIONS 

In the present paper we have studied whether intro¬ 
ducing primordial non-Gaussianities in the density field 
in terms of skewness and kurtosis parameters /nl and 
gNL could have a large impact on the number density 
of nonlinear dark matter haloes at high redshifts. We 
found that values of /nl and gNL much larger than the 
current observational bounds are required in order to 
change the predictions appreciably. Hence, density distri¬ 
butions which are Gaussian modulo skewness and kutosis 
parameters $nl and gNL cannot have a big impact on the 
formation of high redshift super-massive black holes. This 


result contrasts with the large effect which cosmic string 
loops can have m at high redshift. In the case of cos¬ 
mic strings, the distribution is intrinsically non-Gaussian 
and not well described at all by a Gaussian process plus 
skewness and kurtosis. 
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